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The charge current pumped adiabatically through a deformable quantum dot is studied within 
the Green's function approach. Differently from the non-deformable case, the current shows an 
undefined parity with respect to the pumping phase ip. The unconventional current-phase relation, 
analyzed in the weak pumping regime, is due to a dynamical phase shift (f>o caused by the elastic 
deformations of the central region (classical phonons). The role of the quality factor Q of the 
oscillator, the effects induced by a mechanical resonance and the implications for current experiments 
on molecular systems are also discussed. 



INTRODUCTION 



The idea of quantum pumping, i.e. of producing a dc 
current at zero bias voltage by time periodic modula- 
tion of two system parameters, dates back to the work 
of Thoulessi. The pumping is adiabatic when the pa- 
rameters change slowly as compared to all internal time 
scales of the system and the average charge pumped per 
period does not depend on the specific time dependence 
of the parameters. Using the concept of emissivity pro- 
posed by Biittiker et al4, Brouwer— related the charge 
pumped in a period to the derivatives of the instanta- 
neous scattering matrix of the conductor with respect 
to the time-varying parameters. Since then, a general 
framework to compute the pumped charge through a con- 
ductor has been developed for noninteracting and inter- 
acting electrons^. Then the interest in the pumping 
phenomenon has shifted to the experimental^ investiga- 
tions of confined nanostructures, as quantum dots, where 
the realization of the periodic time-dependent potential 
can be achieved by modulating gate voltages applied to 
the structure^. 

On the other hand, the modern miniaturization tech- 
niques have allowed the realization of artificial struc- 
tures of the same length scale of molecular objects so 
that the charge distribution and shape of the system 
may well change during their operation. The process 
of hybridization obtained by joining artificial manmade 
structures with macromolecules makes possible to get 
systems whose electrical response is strongly affected by 
electromechanical tunable interaction. This happens be- 
cause at molecular scale the elastic forces controlling the 
structure of the system are of the same order of mag- 
nitude with respect to the electrostatic forces related to 
the charge distributions. The interplay between electro- 
static and structural degrees of freedom can be experi- 
mentally studied by means of the carbon-based technol- 
ogy allowing a full integration between manmade nanos- 
tructures and molecular conductors (as in particular the 
fullerene Ceo)- In such systems, whose theory has been 
formulated by Gorelik et aL— , the interplay between me- 
chanical and electromagnetic degrees of freedom is re- 
sponsible for a mechanical instability that causes a self- 



oscillating behavior due to a self-consistent bistable po- 
tential. The current induced by this so-called charge 
shuttle phenomenon 9 is proportional to the frequency 
of the oscillating part of the system, while the electric 
charge transported per cycle is an integer multiple of 
the elementary charge — e. The coupling between charge 
transport and vibrational properties of these molecular 
systems has recently been observed by H. Park et aL— in 
an experiment on a single-molecule Cgo-based transistor. 
Similar effects have also been observed in experiments on 
the transport properties of suspended carbon nanotubes 
where the activation of the breathing modes (i.e. radial 
deformations of the nanotube) has been obtained using 
an STM (Scanning tunneling microscope) tipii. These 
systems are nowadays known as nanoelectromechanical 
systems (NEMS) and are very important in clarifying 
the role of the electron-phonon interactions in molecular 
conductors and manmade artificial systems. 
Within the previous scenarios, an interesting question is 
what type of effect would induce the mechanical degrees 
of freedom (oscillations) on quantum pumping proper- 
ties. In particular, one can imagine to realize an adia- 
batic quantum pump with a deformable quantum dot ca- 
pable of changing its configuration under the combined 
effect of electrostatic and elastic forces. This proposal 
can be realized, for instance, by coupling a suspended 
nanotube^ to external leads through two tunable tun- 
nel barriers. Suspended nanotubes of ~ 0.5/im of length 
present a mechanical resonance frequency of ~ 200MHz, 
i.e. of the same order of the adiabatic pumping frequency 
and a quality factor Q ~ 100 — 150 (see for instance Fig. 3 
of Ref.[12]). Thus the mentioned molecular systems are 
good candidates to explore the adiabatic quantum pump- 
ing through a deformable system. 

In this work, by employing a Green's function approach, 
we derive an expression of the adiabatic pumping current 
through a deformable quantum dot subject to a static 
electric field E g generated by a back gate and coupled to 
two external leads kept at the same chemical potential. 
The pumped current in this system is generated by the 
adiabatic variation of the tunneling rates between the 
noninteracting leads and the deformable quantum dot. 
During the pumping cycle the deformable dot is charged 
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with a finite density and reacting to the presence of the 
electric field E g by deforming itself and changing its en- 
ergy. The dynamical effects related to the mechanical 
deformation of the system are considered as adiabatic 
and are described by using a classical formalism. The 
small deformation limit is considered to avoid interesting 
but very complicated non-linear effects^. 
An important point is that the parametric pump we have 
in mind is very different from a shuttle device. Indeed, 
the shuttle consists of a self-oscillating system able to self- 
sustain a mechanical oscillation of a small grain when a 
dc bias is applied to the external leads. In this setup 
the charge is transported on the shuttle and is released 
to the closest lead, the charge on the grain being almost 
constant during the shuttling time. Thus a crucial in- 
gredient is the possibility to reach closely a given lead to 
enhance the tunneling rate which then depends strongly 
on the position of the shuttle in the direction of charge 
transport (longitudinal direction). Alternatively, we are 
proposing here a parametric pumping, the only difference 
from the original work of Thoulessi being the inclusion 
of a classical variable associated to the transverse oscil- 
lations of the central region subjected to a back gate. In 
our pumping setup the tunneling rates are unsensitive to 
the position of the central region in the longitudinal di- 
rection and only depends on the transverse motion of the 
center of mass of the central region. The self-oscillating 
phenomena present in the shuttle dynamics (see e.g. [HI) 
are then absent in our proposal. 

The plan of the paper is the following. In Sec|TT]we in- 
troduce the model Hamiltonian and the relevant param- 
eters. In Sec lII Al the classical dynamics of the variable 
x describing the motion of the center of mass of the cen- 
tral region is studied within the small deformation limit 
(i.e. in linear response). In Sec lII Bl we derive the weak 
pumping formula for the current which includes the ef- 
fects of the additional phase shift 4>d due to the dynamics 
of the ir- variable. In Sec lIIII we present the results of the 
pumped current at zero temperature by discussing the 
role of the mechanical oscillations including the effects of 
the dissipation. Some conclusions are given in Sec lIVI 



II. THE MODEL 

The system under study is shown in Fig[TJ It consists 
of a deformable quantum dot (QD) coupled to two exter- 
nal leads, namely the right lead (RL) and the left lead 
(LL), kept at the same chemical potential fi. The middle 
region is biased by a back gate which generates a time- 
independent electric field E g that exerts an electrostatic 
force —eE e n 7 when the electron charge on the QD region 
is —en (being — e the electron charge). The QD region, 
which can be made of a nanotube or a Cqq molecule, is 
thus deformed under the effect of E g and reaches its equi- 
librium configuration when the elastic force is balanced. 
The electron transport in the QD is generated by the adi- 
abatic parametric pumping consisting in the out-of-phase 



modulation of the two tunnel barrier strength Vi. r (t) at 
the interface with the external leads. 
In order to study the transport through the quantum 
dot we consider as classical the variables associated to 
the mechanical deformation of the QD region, namely 
the position x and the momentum P of the center of 
mass, while the electron transport is treated within the 
quantum theory. Such approximation is justified when 
the center of mass dynamics is much slower compared to 
the electron dynamics (i.e. Vl <C T/ft, being Q, the me- 
chanical frequency of the dot and T the tunneling rate 
controlling the dot-leads dynamics). In the following we 
consider a mechanical frequency f2 of the central region 
comparable to the frequency u of the pump and thus of 
the order of 100 - 300 MHz or less. The Hamiltonian 
describing the system is the following: 

H = ^ tkc\ a c ka + (e - eE g x)d)d (1) 

a=k,L,R 

+ 5^(V a (t)dt Cfca + ff.c.) 

hot 

p 2 Am 2 x 2 

+ 2M + 2 ' 

where the first term describes the Hamiltonian of the free 
electrons in the LL and RL written in term of the cre- 
ation (annihilation) operators c\ a (cka), the second term 
describes the QD energy in terms of the displacement 
x, the third term represents the tunneling between the 
leads and the QD, while the last two terms describe the 
classical degrees of freedom related to the center of mass 
motion of the QD region having mass M and bare reso- 
nance frequency f2. The energy on the QD thus depends 
on the classical state of the oscillator while the oscillator 
is affected by the charge state of the QD. 
Differently from the charge shuttle case, where the tun- 
neling strength V a depends on the longitudinal position 
of the shuttle, in the present analysis such dependence 
is not important and will be neglected in the following 
analysis. 

From the classical Hamilton equations, P = —d x (H), 
x = dp(H) the following equation for the center of mass 
motion^ is obtained: 

Mx + f3x + Am 2 x = eE g {dU). (2) 

where /3x represents a phenomenological dissipative term. 
On the r.h.s. the term eE g {a ! d) is the electrostatic force 
acting on the deformable dot when the charge density 
—e{d^d) is present on it. Under the adiabatic hypothesis, 
the variable x and the ac pumping terms V a (t) can be 
considered as slow varying parameters and thus we can 
compute the retarded Green's function of the quantum 
dot within the average-time approximation^: 

= F r L-r W ( 3 ) 
E — e[x) + i-j 1 

e(x) = e - eE g x, 



where the wide band limit for the self-energy has been 
considered, while T(t) = J2 a 2npV*(t)V a (t), p being the 
density of states of the external leads (pi — p r = p). 
The modulation of the tunnel strengths V a (t) is respon- 
sible for an adiabatic modulation of the linewidth T(t) = 
J2 a T a (t), where T a (t) = T% + T% sin(wf + <p a ). In addi- 
tion, under no bias voltage applied to the external leads 
and within the zero temperature limit (T = 0), the charge 
density —e(cftd) on the quantum dot region can be deter- 
mined by the relation (dU) = tt" 1 Im{G r (E, t)}dE 
whose explicit expression is given by: 

CO" 5 + ^ ctm (w)' (4) 

p, being the chemical potential of the external leads. As 
shown in |4]), the charge density on the dot follows in- 
stantaneously the x-dynamics and thus the equation for 
the mechanical degree of freedom can be solved at each 
time. Once the solution for x = x(t) has been obtained, 
the retarded Green's function of the dot can be calcu- 
lated. 

The instantaneous current pumped in the lead a can be 
computed within the Green's function technique by con- 
sidering the variation of the charge q = —e(cftd) with 
respect to the small and slow harmonic variations of the 
parameters -Xi(t)[i|]: Sq(t) = '}2 i dx i q5Xi(t). By setting 
X\ = e(x(t)) and X 2 — T(t) and by making use of Eq. (|5J), 
the time derivative of q(t) is given by: 

J - -^\G r (E = p,t)\ 2 {T(t)d t e(x) (5) 
+ [p~e(x)]d t T(t)}. 

Noticing that T(t) = J2 a r a(i), the r.h.s. of Eq.©, can 
be written as J2 a ^"(M; t) and thus the expression for the 
current pumped in the lead a within the zero tempera- 
ture limit is recovered^. Under time averaging over one 
pumping period, the first term in Eq.© vanishes and the 
relation F = —I T is verified. The dc current pumped per 
cycle is given by: 

u f 2n/uJ 

ia = ^J Q I a (^t)dt, (6) 

u> being the pumping frequency. In the following we fo- 
cus on the weak pumping limit and confine our analysis 
to the case of small deformation of the central region. 
Within this approximation the equation of motion of the 
x-variable can be linearized and the related solution can 
be expressed in a simple form. 

A. Dynamics of the mechanical degrees of freedom 

In this subsection we treat the dynamics of the in- 
variable within the linear approximation and consider 
the weak pumping limit, <C IV Let us first note that 
(d'd) depends on the instantaneous position x of the de- 
formable quantum dot and such dependence induces a 
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FIG. 1: The deformable quantum dot-based device consist- 
ing of a molecular conductor connected to the external leads 
via two tunnel barrier whose strength can be tuned periodi- 
cally in time by using suitable top gates. In the middle re- 
gion an electric field E g affects the center of mass dynam- 
ics of the QD. Such dynamics takes place in the direction 
perpendicular to the electrons transport. The periodic out- 
of-phase raising/lowering of the barriers strength enables the 
charge/discharge of the central region which reacts to the elec- 
tric field E g due to the electrostatic force. 



non linear term in the dynamics as shown by Eqs.([2]) 
and (HJ. To linearize Eq.([2]), we expand x(t) around the 
equilibrium position x and write x(t) — x + Sub- 
stituting x(t) — x + in Eq.(j4]) we get the following 
weak pumping and small deformation expansion: 

{d U) ~n - ^ Im{G r (E = p)} (7) 

IT 

_J2&l Re {G r (E = p)}, 

* — ' ITT 

a 

where no is defined from Eq.Q by setting £(t) = and 
T%(t) = 0, while Gl{E = p) = [p- e~ + i £ a T%/2]-\ 
e = eo - A2, T%(t) = r^sin(wt + ip a ). Substituting 
Eq.© in Eq.© and by introducing the adimensional 
displacement n(r) = £,/xq and time t — fit we get the 
following linear equation: 

fj + Q- 1 !) + a 2 rj = - A £ K(r)Re{G r (E = p)}, (8) 

a 

where Q represents the quality factor of the oscillator 
(notice that /3 = MQ/Q), A = A/(Affl 2 i ), while a 2 = 
Cl 2 /ft 2 , Cl being the renormalized oscillation frequency: 

fj = ^n 2 + Im{G r (E = p)}\ 2 /(irM). (9) 

The solution of Eq.© can be obtained by the classical 
response function x( T ) which is solution of the equation 
L>x(t) = S(t) where D = + Q^ 1 -^ + a 2 ] is a linear 
differential operator. Under the effect of the forcing term 
F(t), the r.h.s. of the equation ©, the displacement t)(t) 
of the deformable quantum dot can thus be obtained as 
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V(r) = J^dr'xir - r')F(r') 
solution for the x-variable is: 



Therefore the explicit 



x(t) 



^_\X^0 Re{G r iE = fl)}x (10) 



where the quantities \x\ and <j)£> come from the exponen- 
tial representation \\\ exp(i<pjj) of the Fourier transform 
of the response function x( s ) — ( a2 — s 2 — isQ -1 ) -1 , 



lx(s) = 



1 



a 2 - s 2 ) 2 + s 2 Q- 2 



and 



4>d(s) — arctan^-^ — 



From Eq. fTT))) is evident that the dynamical phase <j>jj 
of the classical response function \ plays an important 
role in the quantum pumping. The knowledge of the 
dynamics of the deformable quantum dot allows us to 
derive the retarded Green's function Eq.([3]) and thus the 
pumped current. 



B. Pumped currents within the weak pumping 
case 

In this subsection we derive the expression of the 
charge current adiabatically pumped under the influence 
of the center of mass dynamics of the central region. By 
introducing the retarded Green's function in Eq.Q and 
performing an expansion of the integrand up to second 
order terms in V£ , the current pumped in the left lead i\ 
in units of 



can be written asi£: 

2,77 



-eq 



EM 
2tt 



\G r (E = n)\ 2 Re{G r (E = n)} 



(11) 



x im{G r (E = M )}rf r?[siiu> D ) cos^xrs - r<) 



sm(0D 



r -pi 



s(<f> D )sm(<p)(T l + T5)], 



where the polaronic energy E p — XXxq/2 — X 2 / (2Mfl 2 
has been introduced. The current i l eq is the currem 
pumped when the quantum dot is in its equilibrium po 
sition x and whose expression is: 



L r L l 



sm(ip)Im{G r {E = fi)}Re{G r (E = fx)}.(12) 



Eq. (fTTjl . which is the main result of this work, clearly 
shows that the classical phase <j)r>, coming from the x- 
dynamics, may strongly affect the current-phase relation 
compared to the static situation where an odd behavior 
of the current with respect to the pumping phase tp is 
expected. When A goes to zero, e.g. by switching off 
the central gate [E g — 0), the standard weak pumping 



current is recovered. Furthermore, from detailed analy- 
sis of Eq. one sees that when the pumping frequency 
w becomes very close to the renormalized resonance fre- 
quency f2 the pumped current is significantly modified 
by the additional terms related to the oscillation of the 
quantum dot. Even within the small displacement ap- 
proximation, the amplitude of the second term on the 
r.h.s. of the Eq. (TlTj) proportional to E p Q can become 
comparable to the term proportional to r^T" for suit- 
able values of E g and thus the system may present in- 
teresting oscillation-induced effects on the pumped cur- 
rent. The oscillation- related terms of the current consist 
of three distinct contributions. Two of them are related 
to the pumping phase tp and are proportional to sin0£> 
and cos (/)£>, while the remaining one is independent from 
the pumping phase and can be interpreted as a rectifica- 
tion term. Concerning the terms containing the pump- 
ing phase, the current-phase relation is no more an odd 
function of tp due to the presence of terms of the form 
sin(0£>) cos(v?) and cos(4>d) sm(cp) that can be recognized 
as interference terms. If the oscillator has a very high Q 
(low dissipation), the terms proportional to sin(</>£>) be- 
come relevant only very close to the resonance frequency 
(see Fig.((2])), while the term proportional to sm(ip) is the 
dominant one. Thus for weak dissipation, the odd sym- 
metry of the current-phase relation is preserved almost 
everywhere in the parameters space. On the other hand, 
when the quality factor of the oscillator is decreased, i.e. 
for increasing values of the dissipation, the additional 
terms with phase dependence different from the odd ones 
start to become relevant. Our analysis thus relates the 
dissipation of the center of mass motion to the presence 
of the anomalous terms in the current-phase relation. 
Let us note that the mechanism related to the presence of 
a dynamical phase shift is more general than the context 
we are describing here and can in principle be present 
every time a quantum system is coupled to a subsystem 
with a classical dynamics ( as, for instance, a classical 
RLC circuit). 

Finally, we also stress that the cos(ip) term of the pumped 
current we derived in the present analysis has also been 
obtained in different contexts. For example, in Ref.[l9j 
where the effect of the inelastic scattering within a two 
terminal parametric pump was treated by means of the 
fictitious leads method. In that work a third fictitious 
lead with a time dependent chemical potential was intro- 
duced. In our case, even in the absence of inelastic phe- 
nomena, the deformable quantum dot behaves similarly 
to the fictitious lead of Ref.[19( and thus a cos(tp) term 
appears. Interference-like terms in the current-phase re- 
lation were also obtained in Ref. [2(| where an adiabatic 
quantum pump in the presence of external ac voltages 
was considered. 
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FIG. 2: Behavior of sui(<j!>.d) (dashed line) and cos(0d) 
(dashed-dotted line) as a function of the normalized pumping 
frequency uo. The remaining parameters have been fixed as 
follows: e = 0, E p = 0.015, T l = Tg = 1, Q = 20. 



III. RESULTS 

In the following we fix the zero of the energy at the 
Fermi level (fx — 0) and consider the zero temperature 
limit to study the pumped current (Eq.(fTl|)) as a func- 
tion of the relevant system parameters. All the energies 
are expressed in units of A = 10/zeV which is of the same 
order of magnitude of the static linewidth Tq and the 
polaronic energy E p is assumed as a small quantity com- 
pared to A. Since the polaronic energy can be tuned by 
means of the electric field E g , in the following we assume 
that the displacement of the central region is a fraction 
of 10~ 12 m and thus E p < 10^eV. In starting our anal- 
ysis we notice that the behavior of the pumped current 
depends strongly on the value of the pumping frequency 
since both the response function \ and the dynamical 
phase are functions of u>. This is clearly seen in Fig.® 
where sin(0£>) and cos(</>£>) are shown as a function of 
the pumping frequency to (normalized to the mechani- 
cal frequency fl) while the quality factor is fixed at a 
relatively low value, Q = 20 (high-dissipation). When 
the mechanical dissipation is high, the sui((^d) becomes 
a broad peak close to the renormalized mechanical res- 
onance, while cos(0d) presents a step-like behavior as 
a function of io close to ft. Thus the presence of some 
dissipation mechanism allows for the presence of an even 
function of the pumping phase ip in the current-phase re- 
lation. 

In Figs.Q we analyze the pumped current as a function 
of the pumping frequency to and by fixing the remaining 
parameters as follows: eo = 0.7, E p — 0.015, T l = 1.2, 
T r = 1, T? = 0.3, = 0.5, Q = 150, where we set 
if = ir/2 in the upper panel and tp — in the lower panel. 
The dashed-dotted curves in both the panels represent 
the current computed at E p = 0, i.e. in the absence of 
deformations of the central region, while the dashed line 
represents the current pumped through the system when 
the system is coupled to the classical phonons. Again 
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FIG. 3: Pumped current normalized to —euj/(2n) computed 
as a function of the pumping frequency u. The dashed-dotted 
line corresponds to the current pumped in the absence of pola- 
ronic coupling (i.e. E p — 0), while the dashed curves in both 
the panel represent the current pumped through the system 
in the presence of mechanical oscillations, E p 7^ 0. The up- 
per panel is obtained for ip = 7r/2, while the lower panel for 
ip — 0. The remaining parameters have been fixed as follows: 
eo = 0.7, E p = 0.015, T l = 1.2, F r = 1, Tf = 0.3, = 0.5, 
Q = 150. 



we observe that when the system is driven on resonance 
by the external ac parameters of the pump, the oscilla- 
tions of the dot are responsible for the enhancement of 
the current close to the mechanical frequency Q, while 
for pumping frequencies different from the resonance one 
the system shows a current intensity similar to the static 
case. Furthermore, the specific value of the pumping 
phase f can strongly affect the symmetry of the pumped 
currents as a function of the pumping frequency. Indeed, 
as is shown in the lower panel of Figs. ([3]) for pumping 
frequencies around the resonance frequency f2, the sign 
of the current is reversed. The behavior of the pumped 
current as a function of the pumping frequency u> is q uali- 
tatively consistent with the one measured in Ref . [21( (see 
Fig. 2a of the cited work) where a suspended carbon nan- 
otube with Q — 80 and f2 = 55MHz has been considered 
to analyze the electrical-induced guitar-string-like oscil- 
lation modes. In Fig. ([4J the pumped current is shown 
as a function of the pumping frequency for two differ- 
ent values of the quality factor, namely Q — 75 (dashed- 
dotted line) and Q = 150 (dashed line), while the remain- 
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FIG. 4: Pumped current normalized to — eoj/(27r) computed 
as a function of the pumping frequency lu and by fixing Q = 
150 (dashed curve) and Q = 75 (dashed-dotted curve). The 
remaining parameters have been fixed as follows: ip — n/2, 
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FIG. 5: Pumped current normalized to — eoj/(27r) computed 
as a function of the pumping phase <p. The remaining param- 
eters have been fixed as follows: Q = 150, ui/Q. = 1, eo = 0.7, 
-Bp = 0.015 (dashed line) or E p = (dashed-dotted line), 

v{ = 1.2, rs = i, Tf = 0.15, r? = 0.5. 



ing parameters have been fixed as: ip — n/2, e = 0.7, 

e p = 0.015, r(, = 1.2, = 1, iy = 0.3, = 0.5. As 

shown, higher values of Q (lower dissipation) permit to 
enhance the pumped current close to the mechanical fre- 
quency. The influence of the oscillation-related terms on 



the current-phase relation is shown in Fig. (J5j) where the 
pumped current is reported in the presence and not of the 
deformations. As we see when E p ^0a nonzero pumped 
current appears even for ip = 0, 27r, i.e. a rectification 
term is present. 



IV. CONCLUSIONS 

We studied the adiabatic quantum pumping through 
a deformable quantum dot coupled to two external 
leads in the presence of time varying barriers strength 
and analyzed the oscillations induced effects on the 
pumped current. The out-of-phase adiabatic modulation 
of the coupling to the leads is responsible for the charg- 
ing/discharging of the central deformable region which 
is subject to an electrostatic field E g . The dynamics of 
the center of mass displacement x of the central region 
introduces an additional phase shift (pn which affects the 
pumped current. In particular the deformation induced 
terms modify the current-phase relation by introducing 
even contributions in the pumping phase ip. Such terms 
depend strongly on the quality factor Q and are sup- 
pressed when the mechanical dissipation is low. Further- 
more, the current pumped is enhanced when the pumping 
frequency u> is very close to the mechanical resonance fre- 
quency f2, the enhancement factor being determined by 
the Q factor. Finally, we demonstrated that the pumping 
mechanism can in principle be used to characterize the 
mechanical properties of a molecular resonator allowing 
to study the response function x(i). Thus our proposal 
can be seen as a complementary tool in the investiga- 
tion of the mechanical response of a nanoresonator. On 
the other hand, the deformations of the central region 
permit to enhance the pumped current for values of the 
pumping frequency close the to mechanical one, thus ren- 
dering measurable the pumping current which otherwise 
remains difficult to detect. Experimentally our proposal 
can be realized by modifying the experimental set up pro- 
posed in Ref. [2l[ which is particularly suitable due to the 
relatively low value of the resonance frequency (i.e. 55 
MHz) fulfilling the adiabatic requirement considered in 
our work. 
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